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Abstract 

We consider the issue of an accurate description of the evolution of the non-singlet struc- 
ture function moments M„(Q) near heavy quark threshold. To this aim we propose a simple 
0> , modification of the standard massless MS scheme approach to the next-to-leading QCD 

analysis of DIS data. We apply it to the processing of the modern CCFR'97 data for xF^ 
structure function and extract the value of 



0^ 



O 



as{M-^) « 0.122 ±0.004 
We check also the consistency of light gluino hypothesis with CCFR'97 data. 



^ • 1 Introduction 



XT) 

' An important means of verification of the validity of pQCD (that is, of perturbative QCD 

■ "improved" by the RG summation) is an analysis of deep inelastic scattering (DIS) data. 

To interpret these data within pQCD, one should pay credit to a number of subtle physical 
^ effects: contributions of high twists, nuclear effects, high-order (three-loop) corrections and 

^ ■ the influence of thresholds of heavy particles. All the introduced corrections are roughly of 

^ , the same order of magnitude. 

Qj [ This paper is devoted to the problem of influence of the thresholds of heavy quark (HQ) 

i-S^ ' on the pQCD analysis of DIS data that includes, in particular, the evolution of the strong 

coupling constant as{Q)- Modern estimates performed in [|l|, |2|, ^ have revealed a significant 
^ I role of threshold effects in the as{Q) evolution when the DIS data lie close to the position of 

' "Euclidean-reflected" threshold of heavy particles. The corresponding corrections to Qs(M^) 

can reach several per cent, i.e., they are of the order of the three-loop |^], for the recent 
results see and nuclear effects on as{Mz) ■ 

A common algorithm for the renormalization-group (RG) resummation is based upon 
beta-function P{as) and anomalous dimensions 7(as) calculation and the RG differential 
equations integration performed within the MS renormalization scheme. However, the widespread 
massless MS scheme fails to describe the data near thresholds of heavy particles - b,c quarks 
and, maybe, light superpartners [^. 

An appropriate procedure for the inclusion of threshold effects into the Q ^-dependence 
of OsiQ) in the framework of the massless MS scheme was proposed more than 10 years ago 
0] : transition from the region with a given number of flavors / described by massless 
a.s{Q', f) to the next one with / + 1 ("transition across the threshold") is realized 

here with the use of the so-called "matching relation" for 6ls{Q) 0. The latter may be 
considered as the continuity condition for as{Q) on (every) HQ mass 

a,(Q = M;+i;/) = a,(Q = M/+i;/ + l) (1) 
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that provides an accurate Q;s(Q)-evolution description for Q values not close to the threshold 
region. The condition (0) is used up to the three-loop level; the other version of the matching 
can be found in |p. 

One needs also one more element, the matching procedure for the evolution of the struc- 
ture function moment Mn{Q,m). The corresponding expressions for anomalous dimensions 
7(j)(n;/) are well known in the MS scheme for a fixed / value (see, e.g., ||8|) but until now 
there is no recipe for obtaining a continuous interpolation across the HQ threshold for the 
moment evolution. 

In this paper, we are going to focus just on this aspect of the problem: how does the 
HQ threshold influence the evolution of the DIS structure function? We will examine only 
non-singlet processes of DIS so as to pass over the delicate problem of modification of the 
operator product expansion in DIS through introducing a new scale, the mass of a heavy 
parton (for discussion, see Q). 

To solve the problem we propose a rather simple modification of the massless MS scheme 
to take into account thresholds in analyzing the moments of the DIS non-singlet structure 
function at the two-loop level. 

To simplify the exposition, we shall take advantage of the explicit analytic mass-dependent 
RG-solution derived in |1C] and [11 1 that is expressed directly in terms of as{Q,m) and 



Mn{Q,m) perturbation expansion coefficients. This allows us to avoid the use of RG- 
generators, that are /? and 7-functions. In the next section, we present smooth analytic 
expressions for the evolution of as{Q,m) and Mn{Q,m) at the 2-loop level based on this 
mass-dependent RG formalism. We shall omit all theoretical and technical details (they 



can be found in refs. |10|, |11|, |1^) and write only final results. In Sect. 3, we introduce 
the "spline-approximation" to describe the two-loop level continuous moment evolution, and 
present there another proof of the matching condition (|l|). In Sect. 4, we describe briefiy a 
method of analysis of the DIS data. On its base we carry out the fit of CCFR'97 Collab. 
data, extract the parameter and estimate the contribution of threshold effects. We discuss 
the consistency of MSSM light gluino existence with CCFR'97 data by using the spline-type 
evolution of Mn{Q,m) in Sect. 5. 

Throughout the paper we use the notation: a = as/47r, (a = Qs/47r); indices in brackets 
stand for the loop number, e.g., = Pi-i; instead of the structure function moments 
Mn{Q, m), we consider only its "evolution part" 7W„(Q, m) (i.e., moments of the distribution 
function) 

Mn{Q,m) = Cn{a,Q,m)-Mn{Q,m), (2) 
where Cn are moments of the coefficient function of a certain DIS process (see, e.g., Ref. ||T^ ). 

2 Mass-dependent RG solutions 

In the massless case, the moment two-loop evolution is described by the expression 

M^^HQ) = Mnif^) (3(2^)'"exp{[a-a(2)(Q)] /„} (3) 

with the numerical coefficients 

, _lo{n) f _ Polijn) - (5i-io{n) 

'J'n — „ , In — rO. 1 \^) 



Po = /?(!)(/) = 11 - A/3(i)/ , A/?(i) = 2/3 ; = /3(2)(/) = 102 - A/3(2)/, A/3(2) = 38/3 



In the mass-dependent case, one should use instead of Eq.(^), a bit more comphcated 
expression pH] of the same structure 



DniQ,...) 



exp 



.a(2)(Q,m) 

with the functional coefficients D^, Fn 

r(i)(n,g) 



Dn{Q,m,^) 



^(i)(Q,m,/i) 



Fn{Q,m,fi) 



^(i)(Q,m,/i)r(2)(n,(5) - A(2)(Q,m,/i)r(i)(n,Q) 



and the two-loop running coupling a taken in the form 

1 1 ^(2) (Qj "T-j 

- + ^(1) {Q, m, /i) + - — -— ln[l + a^(i) (Q, m, /i)] 



a^'^^Q , m; a) a 
In the non-singlet case of DIS, 



^m-iQ) = r(i)(n,/) = 7(i)(n)/, / = ln(Q2/^2^ ^ 



(5) 

(6) 
(7) 

(8) 



r(i)fn 

and the HQ-mass-dependent ^(^^ , ^(e=2) Q) appearing in Exp. are just perturbation 
expansion coefficients: 



a{Q,m,iJ.;a)pert = a- a A(^i-){Q,m,^) + a < - ^(2) (Q, m, //) ^ +...; 



MniQ) 



MniQo 



pert 



1 + ar(i)(n, /) + 

+ aH— ^ L+T^2){n,Q,m,f,)\ +... (9) 



satisfying the normalization condition - Ai^£^{Q = /i) = r(^)(n, Q = fi) = 0. 

These coefficients consist of the usual massless part (for / = 3) and HQ-mass dependent 

contributions, e.g., 



1 

/i; a) a 



/3(i)(3)/-A/3(i)^ 

h 



(10) 



with summation over HQ's: /i > 4. Here Ii is the one-loop fermion mass-dependent contri- 
bution, like the polarization operator [15| or the three-gluon vertex loop [^] subtracted at 
Q2 = and C being some subtraction scheme-dependent constant. 

"Massive" RG solutions (^) and (^) possess several remarkable properties: 

• they are built up only of "perturbative bricks", i.e., loop-expansion coefficients ^(f), 
T^£-^{n) (taken just in the form they appear in the perturbative input) and "contain no 
memory" about the intermediate RG entities such as /? and 7 functions; 

• in the massless case with pure logarithmic coefficients, A(^f.-j = r(fc)(n) = 7(fc)(n)/, 
they precisely correspond to the usual massless expressions, like Eq.(^); 

• being used in QCD, they smoothly interpolate across heavy quark threshold between 
massless solutions with different flavors numbers. 



3 Smooth schemes and MS massless schemes 



3.1 Smooth mass-dependent scheme 

We have above considered general formulae to describe the A^„((5)-evolution including the 
threshold effects. It is clear that the mass-dependent MOM schemes automatically provide 
the most natural smooth description of thresholds. To use them in the framework of leading 
order, one needs a mass-dependent expression for Ii presented of Appendix A (for two 
different schemes). So, to perform the one-loop evolution analysis of moments, one should 



substitute Eq. (|2|) or (|2|) in Appendix A into Eq. (0) and then into Eq. (D and Eq. (|), 
and use the approximation 

/ „ \Dn{Q,...) 

We have performed the results of the fit of the "old" CCFR data |19| following this formula 
in l2|. 

However, the MOM scheme meets tremendous calculational difficulties in the next-to- 
leading order of pQCD. Moreover, each of the expansion coefficients Af^^-j, T^^-^ (see Sec. 2) 
becomes gauge-dependent at the two-loop level, which is not convenient. These difficulties 
are absent in the widespread MS scheme. One can go far in loop calculations here (see 
[|13|), but the scheme is not sensitive to the thresholds at all. Below we suggest a practical 
compromise between these different possibilities - the "spline" scheme. This scheme possesses 
both the sensitivity to thresholds and simplicity of the MS procedure. Nevertheless, the MS 
scheme looks like a conventional standard for all DIS calculations now. Therefore, one 
should recalculate the results obtained in other schemes to the MS scheme at an appropriate 
number /. We do not need recalculation for a^^^^^Mz) because the spline and MS (at / = 5) 
schemes evidently coincide at /i = Mz by construction, i.e., a^^^P\Mz) = q^^{Mz; / = 5). 



3.2 MS — vulgate scheme 

Usually, to obtain the evolution law, one calculates numerical expansion coefficients of gen- 
erators P{as), 7n(os)--- in MS scheme and solves the massless RG equations. In the solution, 
with all integration constants being omitted, one arrives at the final procedure which we 
shall name the MS -vulgate scheme. 

The first recipe to include the threshold mass Mf into the framework of the MS evolution 
was formulated in Ref. Q as the "matching condition", Eq. (Q), for the coupling constant. 
Now all measurements on a low scale Q are usually interpreted in terms of the as{Mz) ~ RG 
solution in a certain scheme, with an appropriate matching of different numbers of active 
flavors which evolve from the scale Q to Mz- The matching condition (||) leads to a simple 
rule of including next "active flavors" h into the evolution law 

(0 = /3w(3)^ - A/?(,) • hi ^ A^'p\q, ...) = /?(,)(3)/ - A%Z*; (12) 

i* = J2 [(^{Q' - (^Mo)') HQ/Mhi^)? - (Q - f^) 

h 

Nevertheless, the threshold value M/^^j) does not follow from this procedure and is left un- 
certain. 

Note, the spline-type (in terms of the Z-variable) expression ( p^ ) has an evident anal- 
ogy with the approximation for the mass-dependent MOM scheme formulae for y4(j)((5,/u) 

( r(j)(n, Q, /i)) with the structure A(^i-^{Q,^) ~ f/(j)((5,m) — /(j)(^,m)j, see, e.g., Eq. (|lO|) . 



The approximation being discussed needs an asymptotic form of the mass-dependent calcu- 
lation for the elements I(^i\{z = Q'^/m?), i.e., only logarithmic and constant terms I[i){z) — > 
ln(z) — Cj (see, e.g., Exp. ( p3| ) and ( [27| ) in Appendix A). Based on this form one can construct 
a simple "pure log" ansatz for /(j)((5,m): 

^/Jj^^^^''(Q,M;,(,)) = 0(Q2_M2^.^)ln(QVM2^^^ ; M;,(,) = m^, exp(Q/2). 

This ansatz roughly imitates the "decoupling" property of at Q < M^(j) and 

provides its asymptotic form at Q > It leads to the approximation for j4(j)((5, /i): 

/3(,)(3)/ -A/3(,) J2 [^(Q' - Ml^,)) In (qVM,'(.)) " (Q - ^)] , (13) 

where the threshold position M/j(j) is determined by the scheme dependent constant C(j) . 

A certain value of the threshold M/j(j) in Exp.(|l2|) and another proof of the matching 
condition (||) for the MS scheme can be obtained by using, e.g., the "three-step procedure" 
introduced in Q. Let us review it briefly. 

Well below the threshold, for fj, <C M, one usually uses some effective MS scheme, say 
MSi, that does not take mass of a particle into account. Above the threshold, a new particle 
cannot be ignored, but when Q ^ M, it can approximately be treated as massless within 
some other MS2 scheme. How should the couplings ai(/^) and 02 (Q) in these two MS schemes 
be related? The answer can be obtained by the three step algorithm: 

(i) recalculating from the MSi to MOM scheme at g = <C M to get a*^'^^^(;u); 

(ii) performing the RG evolution of a^^^^ up to g = Q S> M in the MOM scheme; 

(iii) recalculating to the MS2 scheme, including the mass contribution, at q = Q. 

The final result of these successive steps leads to the approximate (due to power corrections) 
equality ~ ruh for the threshold at the two-loop level. We obtain just the same result 
as usually used for the matching condition mentioned above with Mh = rrih- 

Consequently, proceeding in this way, one must modify the perturbative expansion 
coefficients for Ain ,i-e-, r(j)(n,(5) in (|9|), in the same manner as the expansion coefficients 

of the coupling constant = j3il ^ J^^^'^^K For this aim we recall the structure for 

7(2)("'; /) ill the framework of the MS scheme (for details see, e.g., |^], ^ 



7(2) (n; /i + 3) = 7(2) (n; 3) - A/?i • /i • A7(2)(n); (14) 

16 f 3 lln^ + 5n — 3I " 1 

A7(2)(™) = Y |l05i(n) - QS2{n) - - - ^2(^^1)2 | ' ^kin) = ^ 

where the first term ^(2){n; 3) in the r.h.s. of Eq. ( [l4[ ) consists of the usual massless part and 
the parameter h numbers here heavy flavors. It is known (see e.g. |]l^) that 7(2) (^; /) contains 
the terms generated by the evolution of the coupling constant. These terms naturally appear 
in the calculation of two-loop diagrams for ^{2){n; /), they are proportional to the coefficient 
(see the second term in (|l^). Therefore in the MS - expression for r(2) there appears 
a term proportional to the one-loop coefficient A(i-^ : 

r(2)(n, Q; /i + 3) = r(2)(n, Q; 3) - A/3i(/i • /) • A7(2)(n); (15) 



r(2)(n,Q;3) = 7(2) (n, 3)/. 



The /i-dependent part of the ^(i)-term is singled out of Exp.(p^) in the form of A/3i(/i-Z). To 

obtain the continuous coefficient T^^^^\n,Q,h + 3), one should substitute A(^i^ — > A^^^\ 
i.e., {h- 1) — > /* into the second term of the r.h.s. of Eq. (|T5|), according to the recipe (12): 

r(2)(n,Q;/i + 3) ^ rjJ^P^(n, /i + 3) = r(2)(n, Q; 3) - A/^i/* • A7(2)(n). (16) 

Now we can get the complete evolution law by substituting (^) and ( [l6| ) into formulae (^) 
and (|6|), (^) and then into the general formula (|^): 

/ \ A.{Q,M) 

Recent CCFR'97 Collab. experimental data on DIS pO| ] are processed by this method 
(taking also account of the one-loop coefficient function) in the next section. 



4 The QCD fit of the xFs CCFR data 
4.1 Method of QCD Analysis 

In this section, we present the QCD analysis of the CCFR'97 data |2^. They are the most 
precise data on the structure function xF^^XjQ^) . This structure function is pure non- 
singlet and the results of analysis are independent of the assumption on the shape of gluons. 
To analyze the data, the method of reconstruction of the structure functions from their Mellin 



moments is used [21| . This method is based on the Jacobi - polynomial expansion of the 



structure functions. 



Following the method p^, 22], we can write the structure function xF^ in the form: 

A^rr 



' max 



xF,^-^ix,Q') =x"il-xf £ e^'^(x)5]cf (a,/5)Mj^|(Q2), (ig) 

n=0 j=0 

where Q'^^{x) is a set of Jacobi polynomials and c^{a,[3) are coefficients of their power 
expantions: 

n 

@^'Pix) = Y.cf\a,(3)x^. (19) 

j=0 

The quantities Nmax , ol and /3 have to be chosen so as to achieve the fastest convergence 
of the series in the r.h.s. of Eq.(^) and to reconstruct xF3(x,(5^) with the accuracy 
required. Following the results of pl| we have fixed the parameters - a = 0.12 , /? = 2.0 
and Nmax = 12 . These numbers guarantee an accuracy better than 10^^ . 

Finally, we have to parameterize the structure function xF-^ix^Q'^^ at some fixed value 
of Q'^ = Q\ . We choose xFj,{x,Q'^^ in a little bit more general form as compared to |23], 
where the same data have been analyzed within QCD in terms of without thresholds 
effects: 

xF^{x,Ql) = Ax^(\-xf x). (20) 

Here A, B, C, 7 and oq = as(Qo) are free parameters to be determined by the fit. 

To avoid the influence of higher-twist effects, we have used only the experimental points 
in the plane (x,Q2) with 3 < < 200 {GeVjcf and 0.01 <x < 0.75 . The effect of 
target-mass corrections in xFs is taken into account to order M'^/Q'^ p^. 



4.2 Results of Fit and Discussion 



Here we present the results of processing the CCFR'97 data obtained in the framework of 
two different approaches: 

First, we have used the massless /-fixed MS - scheme approach based on the two- loop 
evolution formula @. The corresponding results for ao(/) are collected in the up-parts of 
the of Table 1 for every Q^. 

The second approach is based on mass-dependent evolution Eq.(||), for this formula we 
adapt the "spline" approximation ( p^ . These results for af'' are presented in the down-parts 
of Table 1 for every Qq . 

The results in Table 1 are completed by the results of evolution of ao(/) and af^ to the 
point Q = Mz with appropriate matching of different numbers of active flavors (the last 
column - as{Mz)', see, for comparison with other estimations, review [^]). Repeating the 
fit procedure for different values of Qq = 3, 30, 200 GeV'^ we have obtained the experimental 
dependence of ao(/) on the momentum transfer. In all fits, only statistical errors are taken 
into account. Here we make some comments on the fit results. 

• To demonstrate the sensitivity of ao(/) on the /, the results of the fit are shown for 
different / for each Q^. The largest difference between the values of ao(/) for different 
/ is about 8% in the case of Qq close to kinematical boundaries Qq = 3 and 200 GeV^. 
This difference reduces to 2.5% — 4% variations for as{Mz; f)- There are opposite 
relations for ao(f ) for these two points: ao(3) > ao(4) > ao(5) for Qq = 3 GeV^ and 
ao(3) < ao(4) < ao{5) for Qq = 30, 200 GeV^. Note, the effects mentioned above can 
be described by estimating Aa^ = as{f + 1) — as{f) « das{f)/df : 



dfMnia,f,Q 



exp ) 



daMn{a,f, Qe2 



exp I 



Here, the brackets ((...)) denote the average over experimental values of Ql^p- the 
one-loop level, this expression leads to the simple estimate 

^ ~ -A/3(i)«,((ln(%^))), 

in qualitative agreement with the results in Table 1. 

• The final results for as{Mz) depend on the Qq choice. For the /-fixed scheme this 
dependence is not smaller than 2% (two thousandth of the absolute value) and for 
the spline scheme, it is within only 0.5% of the value of as(M^). Of course, the value of 
af^iQo) lies between the corresponding values of a^^{Qo', / = 4) and a^^{Qo; / = 5), 
that reflects the real distribution of the CCFR'97 experimental points. The suppression 
of the residual dependence on Qq for the spline scheme results gives an additional 
"phenomenological" hint for preferring this scheme over the MS /-fixed version. 

• The results of the fit are rather stable to the mass variations. The 10% change of Mc 
and Mfe yields less than 0.5% change for oq. 

Comparing the spline and MS f-fixed results we arrive at two main conclusions: 

1. The spline scheme is more preferable than the traditional massless scheme, to process 
the experimental data involving thresholds, the values of af'-{Mz,Qo) are focused 
tightly. The average value af^{Mz) = 0.122 is a little bit large than that the CCFR 
CoUab. result for as{Mz) obtained recently in [^] from the Q^-evolution of xF^ and 
F2 for Q^ > 5 GeV"^ and the invariant mass-squared of W"^ > 10 GeV^ 

as{Mz) = 0.119 ± 0.002(exp.) ± 0.004(theor). 



2. The threshold effects reveal themselves as approximately a +1% correction ^to the value 
of as{Mz) ■ The final results for the average values of (5is{Mz) in the fit presented in 
Table 1 look like: 

af\Mz) = 0.122 ±0.001 (t/ieor) ±0.002 (stat) 
af^(Mz;/ = 4) = 0.121 ±0.002 (t/ieor) ± 0.003 (stat). 

Theoretical errors presented here include the uncertainties due to Jacobi polynomial tech- 
nique reconstruction and Qg'deviation of the as{Mz) value in the fit. 



5 CCFR data and the light gluino window 

In Subsec. 4.2, we have obtained a comparatively small effect, one-two thousandth to as{Mz) 
for the threshold contribution at the fit. The order of the effect is determined by the reason 
that only one 6-quark threshold "works" really in the region in question. If Nature would 
provide few thresholds in the experimental interval 3 GeV'^ < < 200 GeV^, then they 
combined influence in a fit becomes significant and the preference of the spline scheme should 
look evident. To demonstrate this here we took an attempt to reconcile the existence of light 
MSSM gluino {g) and the CCFR'97 data. The possibility of the light gluino existence {rrig 
is of the order m^) was intensively discussed few years ago in the context of the discrepancy 
between low energy Og values and the LEP data at the Mz peak. This discrepancy has a 
chance to be resolved by including the light gluino |27]. 

The Majorana gluino leads to large effects in the evolution - = 2 , A/?^2) = 48 

in Eq.(Q) and slows the running both the coupling constant, and the moments A4n{Q)- 
This reinforcement of the contribution of a new g threshold must strongly influence the fit 
parameters. Here we shall suggest the nucleon does not contain light gluino as a parton, and 
gluino reveals itself only in evolution law. It is clear that the standard MS -scheme at a fixed 
/ everywhere is not adequate to the situation. We have performed the fit of CCFR'97 data 
for different values of the gluino mass rrig, the results of the fit for as at nig = 2, 3, 4 GeV 
and rrig = 10 GeV are presented in Table 2. 

Rather a strong growth of aT''^'^\Mz,Qo) value with the decrease of rrig, as well as a 
slight growth of x^i do not provide confidence that the data are consistent with the gluino 
with mass rrig < 10 GeV (see Table 2). Let us consider the gluino mass rrig as a fit parameter 
and "release" it. The best foi' Q^s is reached at m? 

m? = 103+^ Gey 2, a f ^'a\Mz) = 0.138, (21) 

with the asymmetrical stat. errors to the value of m|, (-50)/(±oo) GeV^. The latter 
means that the value of m| obtained in ( |2l| ) is not reliable and the CCFR'97 data "push 
out" the light gluino to the region of a more heavy gluino beyond the fitting region, i.e., 
m| > 200 GeV. Moreover, the "optimum" value of af''^^\Mz) is over 3a higher than the 
world average value for as{Mz) (2^. 

These results are in agreement with the conclusion in that the light gluino is ruled 
out in others inclusive processes. 



^We have found the same estimate, +1%, for the threshold correction to as{Mz) for the processing the 
CCFR date in Q 



6 Conclusion 



We have devised here a new approach for describing the two-loop evolution of the moment 
AdniQ'^) of structure function of lepton-nucleon DIS involving the threshold effects of heavy 
particles. The approach employs an analytic quasi-exact two- loop RG solution |1C, |l^] within 



the framework of mass-dependent remormalization group and on the results of paper [||. 

We adapt here the spline approximation (|l7l) to the above mentioned RG solution and 
obtain as a result the simple modification of the formulae of the standard massless MS evolu- 
tion. For the particular case of a coupling constant evolution this approximation effectively 
leads to the same result as the "matching condition". Finally, this recipe provides a more 
realistic continuous description for the A^n(Q^)-evolution and looks rather simple from a 
practical point of view. We performed the processing of the modern CCFR data to extract 
the value of as{Q) by two different ways: 

(i) the traditional MS-scheme at the fixed numbers of flavors /; 

(ii) the spline scheme with break point at mass M = m; 

The results for the MS spline - scheme processing of the data are the most adequate to 
the situation both for the physical and practical points of view. The threshold contribution 
to the value of as{Mz) consists of about -|-1%; the extracted value of agP^{Mz) is equal 
to 0.122 lb 0.004. The order of the threshold effect is determined by the reason that only 
one b-quark threshold works really in the CCFR'97 experimental region. We examine the 
possibility to reconcile the CCFR data and the MSSM light gluino. The CCFR'97 data 
"push out" the light gluino to the region of more heavy masses. A reconciliation is possible 
for gluinos with mass nig > 14 GeV, but it seems that similar gluinos are less probable due 
to other constraints (see 
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A Appendix 

Here we list the mass-dependent one-loop expressions for Ii{z) and the corresponding /? 
functions in the MOM scheme. First, we write the well-known exact expression for the 
fermion polarization operator { z = /i^/m? ) 

Il^iz) = 2^1 + 4/z (^1 - In (^^1 + z/4 + ^) + ^ - ^ ; (22) 

Iliz ^ oo) = Inz - (c, = ^) + 0(1/ z) ■ Il{z - 0) = ^ + 0{z^) . (23) 
and the /? function 

/3(i) = yC^-^r,^(l-/i,(z)); 
1=1 

where h,{z) = h^z) = ^ - ^^^=|ln (0^+ (24) 

The three-gluon contribution Iq{z) can be expressed in terms of dilogarithms Li2{z), the 
final formula for Iq{z) is too cumbersome and we do not demonstrate it here. As it was 



predicted in the result for the corresponding three-gluon contribution to the /3- function, 
hi{z) = hf{z) -function may be expressed in terms of dilogarithms, as well, 



hUz) = - 



36 



1 



Z ^3/2 ^1 + z/i 
1 1 



In 



^{l + z/3) v/1 + z/4 



In 



1 + 2/4 + Jz/4 + 



1 



+ --^Im [Li2(2;3) - Li2{z^) + -^^2(24) - ^^2(^4)] 



2:3 



1 + ^ 



Z4 



^3/ V2V1 + 2/4 + ^7^, 
The function /if (2;) may be described by the rational approximation [^]: 

{l.2z + l)z 



/if (z) - /if (z) 



(0.15z + l)(0.4z + !)■ 



(25) 



This approximation works well (better than 1% of accuracy) when z > 1, but when z « 
10^2 _ 10-3^ the accuracy is about 10%. One can restore the corresponding approximate 
expression for Iq{z) by the elementary integration of h3{z) 



If{z) = ^•ln(l + z/zi) + (l-^)-ln(l + z/z2) ; 
A = 21/10 ; zi = 20/3 ; Z2 = 5/2 ; 
/f(z^oo) = Inz- {cg^ 2.98) + 0(1/ z) ; 



(26) 
(27) 
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Quark 
content 




2-loops 






ao 


asiMz) 


= 3 GeV'^ 


MS - scheme, / - fixed, Mg = 


uds 


101.2 


0.3674± 0.014 


0.1233± 0.0015 


udsc 


102.9 


0.3545 ±0.0220 


0.1219±0.0020 


udscb 


105.1 


0.3411± 0.0495 


0.1204 ±0.0050 


Sp 


ine, Mc = 1.3 GeV, = 5 GeV 


uds+(c+b) 


103.5 


0.3494±0.0180 


0.1214 ±0.0020 


= 30 GeV'^ 


MS - scheme, / - fixed, Mg = 


uds 


104.5 


0.2156±0.0084 


0.1206 ±0.0030 


udsc 


105.5 


0.2183±0.0080 


0.1214 ±0.0030 


udscb 


106.7 


0.2210±0.0080 


0.1222 ±0.0030 


Sphne, = 1-3 GeV, A4 = 5 GeV 


uds+(c+b) 


105.9 


0.2195± 0.008 


0.1218 ±0.0025 


= 200 GeV'^ 


MS - scheme, / - fixed, Mg = 


uds 


103.3 


0.1632±0.0046 


0.1173 ±0.0030 


udsc 


104.0 


0.1681±0.0050 


0.1197 ±0.0025 


udscb 


104.9 


0.1732±0.0048 


0.1222 ±0.0020 


SpUne, = 1.3 GeV, Mb = 5 GeV 


uds+(c+b) 


104.7 


0.1721±00050 


0.1217 ±0.0020 



Table 1 The results of the NLO QCD fit in the MS-scheme of the New CCFR [gO[ xF^ 
structure function data in a wide kinematical region: 0.0075 <x< 0.75 and 3 GeV^ < < 
200 GeV^ {N exp.p. = 102). The value of the coupling constant is determined for different 
numbers of the flavor (in the up-parts of Table for every Qg) and for few values of the 
momentum transfer = ct(Qo ~ 3, 30, 200 GeV"^). We present in the first column the 
flavor content involved in QCD-evolution using in fit. The results of fit with the matching 
at the thresholds corresponding to = 1.3 GeV and mi, = 5 GeV are presented in the 
down-parts (for every Qq) of the Table. 



New CCFR data and the light gluino window 



Particle 
content 




2-loops 


x' 


"o 




Ql = 2,GeV'\ Mc = 1.3 GeV, = 5 GeV, 


Spline 


M~g = 2 GeV 


uds+(c+b+5) 


111 


0.3102±0.015 


0.1446 


Spline 


M-g = 3 GeV 


uds+(c+b+.9) 


109 


0.3193±0.012 


0.1428 


Spline 


M,j = 1 G(V 


uds+(c+b+ff) 


107 


0.3260±0.019 


0.1415 


Spline 


Mg = 10 GeV 


uds+(c+b+5) 


103 


0.3475±0.018 


0.1372 



Table 2 The results of fit with the matching at the thresholds corresponding to nic = 
1.3 GeV, mi, = 5 GeV and rrig = 2, 3, 4, 10 GeV are presented at the Table 2. The value 
of af''^^\Mz) is calculated with the error about ±0.004. 



